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Abstract 

We study and detail the features of the resonant soliton of the 3x3 operator. The scattering data 
of this operator contains four transmission coefficients, two in each half complex (^-plane, where C is 
the spectral parameter. When the potential matrix has anti-hermitian symmetry {Qa,i3 = ~Q*pa)^ 
each of the two transmission coefficients in the lower half plane become equal to the complex 
conjugates of one of the two in the upper half plane, leaving only two independent transmission 
coefficients, which we take to be those in the upper half plane. The bound state scattering data for 
this operator consists of the zeros of these two transmission coefficients (bound state eigenvalues) 
and a normalization coefficient associated with each eigenvalue. With two transmission coefficients, 
there is a wider variety of possible soliton solutions than in the Zakharov-Shabat (ZS) 2x2 case. 
First, the eigenvalues could belong to only one transmission coefficient, with the other transmission 
coefficient having none. In this case the soliton solution will only exist in one of the two conjugate 
pairs of components of Q lying next to the diagonal, with all other components of Q vanishing. 
These soliton solutions will be equivalent to those in the ZS case, up to scaling factors. If we would 
distribute the eigenvalues among both transmission coefficients so that each transmission coefficient 
would have at least one eigenvalue, the soliton structure becomes more complex. In this case, any 
one eigenvalue will be found to generally make contributions to all non-diagonal components of Q. 
Of particular interest, inside this class is a special class of solitons which exist only when these 
two transmission coefficients have exactly equal eigenvalues. Equality of eigenvalues from different 
transmission coefficients give rise to "resonant solitons". This latter class is the case which we 
treat here. We find that there are two different states for the resonant soliton solution, both of 
which arise from a bifurcation. The bifurcation is shown to arise from the algebraic structure of 
the 3x3 scattering matrix. We also extend the linear dispersion relations (the solution of the in- 
verse scattering problem) to the case when resonant solitons are present. We further show that the 
same result could be obtained by taking the limit of the eigenvalues in each transmission coefficient 
approaching the other. We then study the general soliton solution under anti-hermitian symmetry 
when there is only one eigenvalue in each transmission coefficient, with each unequal to the other. 
We then detail the asymptotics of this solution and then the same when the eigenvalues are exactly 
equal, showing that the latter contains the well known parametric interactions of " up-conversion" 
and " down-conversion" . Lastly, we explain how this equality of eigenvalues in different transmission 
coefficients can be seen to be a nonlinear resonance condition. 



Key Words: Soliton Solutions, Inverse Scattering, Parametric Interactions, Three Wave Inter- 
actions, sl{3) Soliton Solutions, Resonant Solitons. 
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1 Introduction 

The purpose of this paper is to detail the features of the resonant sohton solutions of the nonde- 
generate 3x3 eigenvalue problem 
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on the interval —00 < x < +00. We assume Ji > J2 > J3 and that Tr(J) = 0. Q{x) is a potential 
matrix which vanishes like Q{x — > ±00) = o{l/x) for large x. The matrix V{x) is a 3 x 3 solution 
matrix which contains Jost solutions as its columns. This eigenvalue problem provides the Inverse 
Scattering Transform (1ST) for the three- wave resonant interaction (3WRI) |47l I24|, H8 l I25|. Here, 
we generally assume no symmetry for the matrix Q and will take the six components of Q to be 
generally independent and uniquely different. 

The Lax pair for the 3WRI nonlinear system was first presented in 1973 [47] along with the 
soliton solutions for the explosive and decay cases. Here the components of Q are the envelopes for 
interacting waves. Here for the first time, it was demonstrated that in the explosive case, there was 
a nonlinear instability whereby a singularity with infinite amplitude would eventually develop in the 
solution. In 1975, approximately at the same time, two manuscripts were submitted [241 148j giving 
more details on this system and its solutions. The work by Kaup pi] presented a derivation of the 
1ST for this system, giving a set of Marchenko equations and described the solutions from the point 
of view of the Marchenko equations. In the work by Zakharov and Manakov [l8], they constructed 
various soliton solutions and described the interactions from the point of view of these soliton 
solutions and the quasi-classical approximation. Each work made use of separable initial data for 
constructing the general scattering matrix. The final solution and the entire interaction could then 
be described in terms of the initial soliton and radiation (continuous spectra) content of each wave. 
In 1979, Kaup, Reiman and Bers [25\ publish a review of the 3WRI wherein more detailed aspects 
of these solutions were given along with comparisons and validations from numerical solutions. In 
a related work [35], Reiman extended these results to the case where the medium could contain 
spatial inhomogeneities and studied how they differed from the homogeneous solutions and how 
they modified those results. Interestingly, Reiman showed that the inhomogeneous case of the 
3WRI was also an integrable system solvable by the same 1ST of the homogeneous 3WRI. Higher 
order 3WRI soliton solutions, particularly those corresponding to multiple zeros in the transmission 
coefficients, have recently been given in Ref. [39]. A perturbation treatment of the decay instability 
was briefly treated in [25]. More recently, a full treatment of the perturbations of (II. ID and the 
universal covering set for the squared eigenfunctions and their adjoints has been given in Ref. j30j . 
including the closure and the orthogonality relations. The general nature of the soliton solutions of 
the n X n generalization of (jl.ip has been described in [50) . These soliton solutions have also been 
classified in Ref. [20| according to the group properties of the eigenvalue problem, while a recent 
overview of the n x n 1ST and its soliton structure is found in Ref. 

For the simple 2x2 case [45^ [T], one has only one type of soliton, generally called an "NLS 
(Nonlinear Schrodinger)" soliton. Due to the group properties of the Zakharov-Shabat (ZS) eigen- 
value problem, the NLS soliton is classified as a sl{2) soliton [20]. The soliton solutions that arise 
in the 3x3 case are classified as sl{3) solitons. The general soliton solutions for three of the sl{3) 
symmetries have been given in [A8\ [2^ 125]. 

In the 3WRI, the upper triangular components of Q correspond to three different waves which 
can interact resonantly. (The lower triangular components of Q are it the complex conjugates of 
those in the upper triangle.) There are two general types of s/(3) soliton solutions. First, depending 
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on the symmetry, any one of the three waves could contain one or several NLS solitons, provided all 
the solitons were found in only one wave. This is a trivial example of the sl{3) solitons wherein one 
has only one nonzero wave, not three. Thus no interactions occur. These sl{3) solitons are identical 
to sl{2) iV-soliton solutions. This is expected since sl{2) is a subgroup of sl{3). Thus there would 
be three different representations of the s/(2) subgoup, one for each wave. The second general type 
is the nontrivial case, wherein all three waves are nonzero. For these s/(3) solitons, they interact 
with their energy being exchanged back and forth between the three waves in a regular fashion. 
However the asymptotics of these solutions show that in general, no solitons are exchanged between 
the three waves. 

In the general 2x2 eigenvalue problem where the two off-diagonal components of Q are in- 
dependent and unrelated, there are two transmission coefficients, one in each complex half-plane, 
the zeros of which, along with a normalization coefficient associated each zero, determine the sl{2) 
soliton structure. In the case of the 3x3 eigenvalue problem, one has a total of four transmis- 
sion coefficients, two in each complex half-plane. The zeros of these, along with a normalization 
coefficient for each, determine the s/(3) soliton structure. Amongst these structures lies one that 
is uniquely different from sl{2) solitons, which is called a "resonant soliton" |30l I31| . A resonant 
soliton is one where the two transmission coefficients in each half-plane have a common zero which 
gives a bifurcation in the general sl{3) soliton case mentioned above. It is unique among the sl{3) 
soliton solutions in that not only is energy exchanged between the waves, but also entire solitons 
can be exchanged. This feature of the 3WRI, whereby one can transfer packets of energy from a 
wave of one frequency to another wave of a different frequency, and/or back again, has long been 
suggested as a basis for possible designs of optical logical elements. 

Recently this feature has generated further interest for the same purposes in a different setting. 
Here one utilizes a nonzero background wherein one takes a combination of bright and dark solitons. 
Degasperis, et al [H [T^ [T3] have found exact solutions of this 3WRI when one wave is taken to 
be on an asymptotically constant nonzero background, with the other two being localized pulses. 
Solutions exist where the three interact and move with a common velocity (simultons). Such 
solutions were found to be stable when their velocity was greater than a critical value. Interactions 
of simultons with various localized pulses were found to give rise to the excitation (decay) of stable 
(unstable) simultons by means of the absorption (emission) of the energy carried by a localized 
pulse. The speed of these solitons could be continuously varied by means of adjusting the energies 
of the two bright pulses. 

Here we shall only treat the bright soliton case. The major purpose of this paper is to detail 
those features of the scattering data which allow this soliton resonance phenomenon to occur. 
Secondary will be to describe how these features will display themselves in the potential matrix, Q. 

As is known, physical applications of the 1ST typically require the potential Q to have certain 
symmetries. In the 2x2 case on the infinite interval, the symmetry for most physical applications 
are found in one of the two cases Q21 = r = ±q* = ±Q|2 where * indicates the complex conjugate. 
If Q{x ±00) — 7> 0, then solitons can only occur in r = —q* case [45^ [l]. If one allows for an 
asymptotically constant nonzero background, then in the r = +q* case one can have "dark solitons" 
|46j . If the solution is to exist only on a semi-infinite or finite interval, virtual solitons [41 ^ I27j for 
either symmetry can be found. 

For the 3x3 case, one has several other possible symmetries, three of which are found in the 
3WRI [48t Y2M, [25] . Here we shall only consider the general case of no symmetry or the symmetry 
Qoi,i3 = ~Q*^ Q,) which we shall call the anti-hermitian symmetry. This latter symmetry in the sl{3) 
case corresponds to the soliton decay case of the 3WRI and in the sl{2) case, to r = —q*. 

In Section 2, we will briefly outline the results from the direct and the inverse scattering problems 
for Eq. (jl.ip . For the case of no symmetry, we give the linear dispersion relations (LDRs) from 
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which one defines the scattering data. For the anti-hermitian symmetry, the bound state scattering 
data consists of the zeros of two different transmission coefficients, with a normahzation coefficient 
associated with each zero. Prom this structure, as discussed earher, there are two general types 
of sl{3) soliton solutions. First, zeros could be in only one transmission coefficient, with the other 
transmission coefficient having no zeros. In this case, one has only an sl{2) A^-soliton solution which 
would only exist in one of the two conjugate pairs of components in Q adjacent to the diagonal. 
For the second type, the zeros would be distributed between both transmission coefficients. This 
soliton solution would make contributions to all three conjugate pairs of Q. To specify these s/(3) 
soliton solutions, we use the notation (M, A^)-soliton solution, where M will be the number of zeros 
of the first transmission coefficient and will be the number of zeros in the second transmission 
coefficient. We also give and discuss the general s/(3) (1, l)-soliton solution. When these two 
eigenvalues are equal, we have a resonant soliton. 

In Section 3 we discuss the general s/(3) soliton solution when an arbitrary number of resonant 
solitons are present. We find that there are two distinctly different structures for a resonant soliton, 
due to a bifurcation in the scattering matrix. These structures are obtained for each branch of the 
bifurcation, for the case of compact support. Also the LDRs are extended to include resonant 
solitons. We then briefly discuss how to handle the non-compact support case and how it will give 
the same structures for resonant solitons. In Section 4, under anti-hermitian symmetry, we obtain 
the resonant sl{3) (1, l)-soliton solution from the extended LDRs. Within this solution, we flnd 
the third representation of the sl{2) subgroup and delineate its scattering data. In Section 5, we 
briefly discuss how the resonate soliton can also be obtained by taking the limit of the two zeros 
approaching each other in the general sl{3) (1, l)-soliton solution. The same can be also be done to 
obtain the extended LDRs. A summary of the results and conclusions is given in the last section. 

2 The Direct and Inverse Scattering Problems 

Here we highlight the direct and inverse scattering problems for the general 3x3 operator given in 
[30] , which followed the method used in Refs. [171 [Ml SSI |25] . This method allows one to obtain a 
linearly independent set of the LDRs, from which one could reconstruct the potentials. Later, we 
shall employ these LDRs to obtain the general (l,l)-soliton solution for the 3x3 operator. 

2.1 The Direct Scattering Problem 

In the direct scattering problem, one addresses the solutions of the eigenvalue problem, what their 
analytical properties are, the adjoint solutions and their properties, what is the scattering matrix 
and its properties, what are the features of the bound states, if any, and what are the fundamental 
analytical solutions and their adjoints, etc. Each of these topics shall only be defined and discussed 
to the extent necessary for the construction of soliton solutions. Further details will be found in 
the above references, gZl [Ml SSI [25l [30] . 

The relevant matrix form of the linear eigenvalue problem associated with the 3WRI is given by 
(jl.ip . We assume no symmetry on Q and take the six components of Q to be independent (all the 
diagonal components are zero). We shall take the diagonal elements of J to be real and to satisfy 
Ji > J2 > J3. Qix) is a potential matrix with zero diagonal entries and vanishing as \x\ 00. 
The matrix V{x) is a 3 x 3 solution matrix which contains the Jost solutions as its columns. For 
(jl.ip . we shall define the Jost solutions by their asymptotics as x ^ ±00. For ^ real, there are two 
standard sets which are 

$(x ^ -00) ^ e''^-^^, ^{x ^ +00) = e^^^^ . (2.1) 
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Since these two matrices are linearly dependent on the other, we have 



$ = ^' • 5, ^ = ^-R, 



(2.2) 



where due to the Wronskian relation, 



detS" = 1, detR = 1, 



while the inverse relation is 



R 



S33S22 — S32S23 »S'325i3 

5'3l523 — 533521 5335ii 

00 O C CO 

»-'32>-'21 — '-'31'-'22 '-'31'-'12 



5335l2 
53l5i3 
5325ii 



523^12 — 522^13 
52l5i3 — 5235ii 
5225ll — 52l5l2 



(2.3) 



(2.4) 



We will use it superscripts to indicate the regions of analyticity (+ for the UHP and — for the 
LHP) where UHP stand for the upper half complex (^-plane and LHP for the lower half. Subscripts 
will used to indicate the components of a matrix quantity. The individual Jost solutions are the 
columns in the Jost solution matrix. For this system, we have 



[V'^",'02,V'3 



(2.5) 



where those components without a it superscript in general only exist on the real (^-axis. (Strictly 
speaking, it is not the Jost solution which is analytic in C, uniformly for all x, but the various 
columns in matrix products such as $ • e~*^"^^. With this understood, we shall refer to the ap- 
propriate Jost solutions as being analytic in C if this product is so analytic.) How to generally 
determine the analytical properties of the Jost solutions are detailed in the above references and 
textbooks such as Ref. [38} [50]. 

The analytical properties of the scattering matrices S follow from the above and are 



(2.6) 



In each region of analyticity, there are three linearly independent solutions of (11. ip . For the real 
axis, those could be either $ or ^, or a suitable mixture. For each half-plane, we already have two 
of these Jost functions. The third can be construct from a linear combination of the Jost functions 
as detailed by Shabat [38]. Such solutions were constructed in [2l]. For inversion about x = +00, 
one can take our third independent function to be the meromorphic functions |30] 
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We define the set of solutions, /u^ and ^ , as 



(2.7) 



(2. 



From these solutions, one may construct the "fundamental meromorphic solutions" (FMS) 
which will be designated by O. They are defined by 



{2S 



0"*" provides us with a set of three linearly independent, meromorphic solutions in the UHP and 
existing on the real C-axis, while Q~ provides us with another similar set on the real ("-axis and in 
the LHP. 
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The remaining part of the direct scattering problem is to detail the asymptotics of the Jost 
solutions as one approaches any essential singularity on the boundary of the region of analyticity. 
This provides a means for obtaining the potentials, given the Jost solutions. There is only one 
essential singularity at | C |= oo in this problem. Taking to be real, then $ and ^ have a common 
asymptotic expansion which is 



$ , ^ = (I3 + /Q + /C^ + ...)■ e-^^-^^ as Id 
In terms of 0^, this leads to 

G±(|C|^oo) = l3 + i^«/C + 0(l/C'), 



00 . 



(2.10) 
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in the appropriate half-plane. One finds that when no pairs of the components of J are equal, then 
^(1) 

can be given by J] = Q, whose solution is 
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where the X's are generally integrals over quadratic products of the components of the potential 
matrix, Q. 

2.2 The Inverse Scattering Problem 

On the real axis, which is the boundary between the two regions of analyticity, it follows that 
and 0~ will be related linearly. Thus one may construct what are called "linear dispersion 
relations" (LDRs), whereby one obtains a set of nonhomogeneous, linear, algebro-singular integral 
equations relating 0+ in the UHP (and on the real axis) to Q~ in the LHP (and on the real axis), 
and vice versa. As shown in |30j . one can reduce these six equations to a set of three linearly 
independent equations, from which one can define the scattering data and which also provides a 
solution for the inverse scattering problem, if it exists. (As is well known, sometimes potential 
components in Q have singularities which are physical, as in the explosive case of the 3WRI [47J.) 
An alternate derivation of the LDRs has been given in Ref. [31], whereby one directly obtains a 
set of 3 linearly independent LDRs for the 3x3 case of (jl.ip . which is obviously extendible to the 
problem of obtaining a set of n linearly independent LDRs for the general n x n case of (11. ip . 

Taking the results from ^30j, we have the four LDRs (of which only three are linearly independent 
on the real axis) 
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and for 9(0 < 0, 
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c) 



^''7^0^(C3\fc)i^23(C3^3..) 



(2.18) 
(2.19) 



33, fc 



fc=l (^33,fc 

As to the notation and quantities given above, TZ indicates that the path of the integral is along the 
real axis, N^^ (-^33) is the number of zeros of Rfi{C) ('S'33(C)) iii the UHP (assumed finite), Ctik 

C33 fc ) is the A;th zero of Rii{C) ('S'33(C)) i^i the UHP, N^^ (-^33) is the number of zeros of ^^^^(C) 



(i?33(C)) in the LHP (assumed finite), Cii ^ (C33 fc) is the A;th zero of 5'-^-|^(C) (-R33(C))i k l*^' 



-^33, fc ' 

"^n fc' -^33'fcj i^ value of dR'l^{C,) / dC, [dS^^{C)/dC, dSi^{C)/dC, dR^.^{C)/dC) at its fc th zero, and 
-Ri2,fc) etc., in the case of compact support, are just the values of R12, etc. at the appropriate zeros 
(and in the case of non-compact support, would just be some coefficients). Finally, 



EpqiC) = exp [iCiJp - Jq)x] . 



(2.20) 



We observe that the integrands of the integrals in (|2.16p and (|2.17p are equal and oppose in sign 
while the discrete contributions in each are identical. From their difference, taking the limit of the 
imaginary part of C vanishing, we have the relationship 



-R32 
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©J - ©2" = ©3^ ^32(0 --^QJ EuiO , 9(C) = , 
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33 



(2.21) 



which according to (12. 7p . is an identity. Whence ()2.16p and ()2.17p . for ^J(C) = 0, are linearly 
dependent. However, since the other two equations require the values of 0^ at the various poles, 
we need to retain both of these in the set of LDRs. 

This total set of nonhomogeneous, linear, algebro-singular integral equations is a minimal set 
of LDRs, from which one may reconstruct the Jost solutions, given the scattering data. 
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2.3 The scattering data 



Once we have the LDRs in the above form, then it becomes possible to define the scattering data 
for this problem. In order to solve these equations, we first must specify the quantities: 

• the reflection coefficients aji = {j = 2,3), cr^s = (j = 1,2), pi2 = — tt and 

^33 -^11 

P32 = — on the real axis, 

^33 

• the zeros of Ri^iC) in the UHP {C+^ j^; k = l,2,...,iV+) and the values of Ci2,k = at 

' ' ^11, fc 

each such zero, 

• the zeros of S'^(C) in the UHP (^3^ ^; k = 1, 2, A^'3'3) and the values of C23,k = ir^ at each 

' ' "^33, fc 

such zero, 

p 

• the zeros of R^^{C) in the LHP (Casfe; ^ = li 2, A^g^) and the values of C32^k = j^-i'' at 

' ' ^33, fc 

each such zero, 

S 

• the zeros of «S']j^(C) in the LHP ((^7, ■ ; k = 1, 2, A^j^) and the values of C21 k = at each 

' ' '^ii.fc 

such zero. 

Note that we have six reflection coefficients but only four sets of eigenvalues and normalization 
coefficients. Under compact support, the four normalization coefficients would be the residues of 
four of the reflection coefficients at the appropriate eigenvalues. Thus there are two reflection 
coefficients, (T13 and CJ31, which will not be associated with any normalization coefficients. In the 
linear limit, each reflection coefficient can be matched to the Fourier transform of one of the six 
components of Q. With only four normalization coefficients, we may specify the position and phase 
of solitons in only four of the potential components. If solitons are to occur in the remaining two 
components, then their positions and phases cannot be given independently. Instead the positions 
and phases of solitons in these components must be determined by the positions and phases of the 
other four solitons. It is this point that we want to further understand. We shall return to this 
point later. 

Under the anti-hermitian symmetry, the adjoint solutions of (jl.ip are linear combinations of the 
hermitian conjugate of the Jost solutions, V . As a consequence of this symmetry, Ra^p{C,) = S*^ ^(C*) 
which, on the real (^-axis, leads to a2i = PI2 and P32 = (723. For the bound state data, N^^ = N^^, 

^33 = ^33' Cn,k = Cntfc' C33,k = Ct3*k^ C2i,k = C*2,fc and C32,k = C^23,fe- Thus each quantity in 
the LHP is the Hermitian conjugate of another quantity in the UHP. This symmetry gives rise to 
nonsingular s/(3) soliton structures which most closely match that of the familiar sl{2) solitons 
where r = —q*. In the following subsection, we shall obtain the (1, l)-soliton solution under this 
symmetry, which is the simplest nontrivial s/(3)-soliton solution. Using that, we will describe the 
resonant nature of this soliton solution. 



2.4 The (l,l)-soliton solution 

The solution method for the (M, A^)-soliton solution when there is no symmetry, is given in Ap- 
pendix A. We use those results to obtain the (l,l)-soliton solution for the anti-hermitian symmetry, 
provided that the two eigenvalues are distinct. To that end, take A^^ = N^^ = N^r^ = N^^ = 1. 
Since there is only one zero of each transmission coefficient, we can simplify the notation. We elim- 
inate the triple subscripts on the eigenvalues, replacing them with only one subscript, replacing 
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Cfi I with C^, with (^^ , etc. We also drop the last subscript on the C"s since there is only one 
eigenvalue arising from each transmission coefficient. 

Turning to ()2.1ip and ()2.12p . we may solve for the components of Q. First we impose the 
symmetry condition Qa/s = —Q*^^ which translates into the symmetry Rai3{C) = 'S'^a(C*) ^or the 
scattering matrices. There are four fundamental quantities in the LDRs, on which the solutions 
depend. These are given below and due to the assumed symmetry, can be reduced to only two 
x-dependent complex quantities. 

7^^2= ^Ei2,+i= e^^i2 gifi+(Ji-J2)(x-xi2) ^2.22) 



2r?i 
C21 

TZ32= £^32,-3 = e-^^23e-<3"(-^2--^3)(x-r.23) (2.24) 



521 = fj- £^21,-1 = (2.23) 



523 = §f i^23,+3 = ^32 (2-25) 



where r]i is the imaginary part of Ct^ % is the imaginary part of ^3^, 712 and 723 are real phases 
and X12 and 3:23 are real spatial positions. Then from the results at the end of Appendix A and 
equations ()2.1ip and (|2.12p . we have 



2ir?i^i2 A , (3" - C+ 



gi2 = (Ji - J2) ^^^^ [I + |7^32|'J , (2.26) 

Qi3 = (^i-^3)%5^, (2.27) 



^(C3^-Ci+) ' 

2^%^ f, , C3+ - cr ,^ ,2 



Q23 = -(J2 - J3) ( 1 + ] , (2.28) 



where 



D=(l + |7^l2p) (1 + |7e32p) + , ^"^^^2 1^12!" 1^321' , (2.29) 

^ ^ IC3+-C1+I 

and the other three components of Q follow from the assumed symmetry. This solution was 
originally given in Refs. \A7\ although in a different form. 

Contained in these solutions are the two simple (1,0)- and (0, l)-soliton solutions. For the 
(1, 0)-soliton solution, we would have 5*3^3 without a zero and the residue 7^32 would be absent from 
(j2.26p - (j2.29p . In this case only Q12 would be nonzero and it would have the form 

= (Ji - J2) ^^T?^ ' (2.30) 

1 + |7Ci2| 

which is proportional to a simple sl{2) soliton located at x = xi2- Similarly for the (0, l)-soliton 
solution, we would have R^^ without a zero and the residue 71x2 would be absent from (j2.26p - (|2.29p . 
Now only Q23 would be nonzero and it would be given by 

Q23 = -(J2-J3)-^j^, (2.31) 
1 + I /V32 1 

which is also proportional to a simple sl{2) soliton located at 2; = X23. 

Note that there are only two positions and two phases in (I2.22p - (I2.25P while there are three 
waves. Thus one wave can be considered to be determined or driven by the other two. If we take 
Q13 to be driven, then its phase and position is determined by Q12 and Q23- As to the variety of 
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the solutions that one could obtain here, one finds both analytical and numerical solutions in the 
literature (see, for instance, [8], [32], [13], [33], [l2], [T5]). 

Before leaving this solution, there are features of it that we should point out, particularly its 
asymptotics. To do this, we shall take iCa^ — Ci'l to be generally nonzero and on the order of the 
imaginary parts or even larger, and the differences Ji — J2 and J2 — J3 to be roughly the same. Let us 
translate our x coordinate so that one of the two positions will be at zero. Then the general nature 
of this solution can be obtained by letting the other position vary from —00 to +00. The phases 
are relatively unimportant in determining this nature so they will be ignored. Taking X12 = 0, 
then for X23 large, either positive or negative, Q12 is an NLS soliton, generally localized around 
X = 0. Meanwhile Q23 is also an NLS soliton, but generally localized around x = X23, while Q13 
is generally exponentially small. As |j;23| decreases (the two NLS solitons approaching each other) 
and approaches zero, the two NLS solitons interact and Q13 grows. However its growth is limited 
by two main factors. First, if the real parts, or even the imaginary parts, of the eigenvalues are 
widely divergent, then the maximum growth in Q13 is obviously limited by the |(^^ — Ci"! term in the 
denominator. Second, the mathematical structure of the solution for Q13 is such that the amplitude 
of Qi3 will generally be bounded by about a third of that of an NLS soliton with a similar width. 
In terms of the 3WRI, the variation of X23 from —00 to +00 models two NLS solitons (in different 
waves) coming together, colliding, interacting, and then during the interaction some amount of Qis 
is produced. As the two NLS solitons pass through each other and then recede, what energy was 
in Q13 is returned to the two NLS solitons, which have suffered no consequences for this collision 
except for shifts in positions and phases. This regime of the (1, l)-soliton is a near-resonant regime 
wherein only a brief temporary conversion occurs. 

Let us now take up the singular dependence in (j2.26p - (|2.29p . which depends on the difference 
of the two eigenvalues, which up until now, we have taken to be non-zero. If one takes TZi2 and 
7^32 to be fixed and take the limit where this difference vanishes, then all components of Q vanish. 
However this ignores the dependence on x. As x increases, both 7^12 and 7^32 will smoothly pass 
toward zero. One concludes then that as the difference in the eigenvalues approach zero, for fixed 
C12 and C32, the dominate change to the soliton structure will be that the entire structure will 
simply shift to larger values of x. Alternately, one could achieve the same effect by allowing the 
amplitudes of C12 and C32 to become smaller and smaller, and at such a rate that the general 
position of the soliton structure would not shift to larger x values. However there is a problem if 
we allow C12 and C32 to vanish: we lose the soliton phase and position information. So there is 
more involved here than simply a scaling. 

To understand this better, for the sake of argument, let us consider the case of compact support 
where all components of the scattering matrices can be extended into the entire complex plane. 
Then from (j2.4p . we have 

^33 ~ ^11 ^22 — R12R21 ■ (2.32) 

Now, R'^^{(^) = and S^^{C^) = 0. So if C3'^ Ci^ then the left hand side vanishes as well as the 
first term on the right. Thus it follows that if Cs^ = Ci", then the product R12R21 must also vanish 
(whenever the potential is on compact support). Now, compact support is crucial to this argument 
and pure soliton solutions are never on compact support. Nevertheless this still indicates that care 
will have to be used whenever these two eigenvalues approach each other. We will resolve this in 
more detail in the next section. 
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3 The resonant soliton case 



With the above in mind, let us study this particular singular case where these two transmission 
coefficients have exactly equal zeros. As one may easily verify, when there is a zero of Rii{C) iii 
the UHP, ^, which exactly matches one of the zeros of 'S'^(C) in the UHP, Css^; oi' when there 
is a zero of R^^iC) in the LHP, ^, which exactly matches one of the zeros of <S'{'^(C) in the LHP, 
£, then one will find a singularity in the LDRs, (|2.13p - (j2.18p . For example, let us take Cn k to 
be exactly equal to C33 £ for some value of k and i. Then to solve these LDRs, we need the value of 
^tiCtse) *° insert into (I2.19p . This value will follow from (I2.16p . However, when we evaluate this 
equation at C = Cas^j we see that there will be a term of the form 



Rl2,k 

+1 



which has a denominator that, if -Ri2,fc0i'(Cii fc) nonzero, becomes singular if Ctik exactly 
equal to Css^- Thus if a solution is to exist, Ri2,k^i {Cti k) ™ust vanish. That is essentially what 
happens, which we shall now detail. For simplicity, we shall assume compact support so that R 
and S*, and the Jost solutions $e~*^'^^' and ^e~*^'^^' are entire functions. Later on in this section, 
we shall briefly discuss how to handle the noncompact support case. 

First, we will collect the consequences of having such paired eigenvalues and introduce our 
notation. Since R22 and 5*22 are not transmission coefficients, we have the subscript "2" free to use. 
So let us designate this common zero instead by C,2k where the range of k will be over all such zeros 
in the UHP. We shall assume this zero is simple in both transmission coefficients. Consider the 
(1, 1) component of (|2.4p at = Ctk- Then it follows that the product S'235'32 must also have a zero 
at C = ctk- -'■^ order to avoid double zeros in this product, we shall also assume that 8228^^ ^ Rfl 
at C2k- compact support, either 523 or ^32 must contain this simple zero. If we take ^23 

to have this zero, it follows from (12. 4p that R21 must also have this same zero. On the other hand, 
if we take 5*32 to have this zero, then it similarly follows that R12 must also have the same zero. 
Thus we find that we have a bifurcation with two different possible options. 

Similarly in the LHP, assuming no symmetry, when C = C2k ^ simple common zero of both 
i?33 and S^i in the LHP, then it follows that 512 5*21 must also have the same simple common zero. 
Taking S21 to be non-zero at C^^, we find that S12 = = i?32 at ^2~fc- If we take 5i2 to be non-zero, 
then S'21 = = i?23 at (2,k- 

These conditions follow directly from the matrix structure of the sl{3) group and compact 
support. All that has to happen is for a zero of Rfi{C) to match some zero of 5';^(C) in the UHP, 
with both zeros being simple. Then the other zeros follow, and similarly in the LHP. 



3.1 Bound state structure - UHP 

In the following, we shall assume no symmetry on the matrix Q, which we shall take to be on 
compact support. Let us turn to the structure of the Jost solutions when Rf^ and have 
a common zero. In the following, we shall freely make use of ()2.4p and the two orthogonality 
relations (in tensor notation) 

Ra,lSl,(3 + Ra,2S2,l3 + Ra,3S3^(S = 6'^ = Sa,lRl,l3 + Sa,2R2,l3 + Sa,3R3,l3 , (3.1) 

where 5? is the Kronecker delta. 



11 



First, from our definitions of /u^ in (j2.8p . we have that /i^^ , Rf^ — xt ^ -^^3/^2 ~ X2 ^^"^ /^3 ' 
are entire functions of where X2 ^^^^ taken to be defined by the above relations. In general 
jjL^ is a meromorphic function while is an analytic function [30| . Consider the T-matrix and its 
inverse given in Ref. |30j 

^l+ = ^L-■T, fi- = fi+-T-\ (3.2) 

where 





- 1 


R12 


Sl3 




- 1 


5*12 


^13 




5*11 


SiiRii 






Rn 


^11^33 




T = 


5*21 


^ 1 R12S21 
SiiRii 


R23 




R21 

•^33 


-R325'23 

Q+ P- 
'-'33-"'33 


5*23 
•^33 




R31 


S32 


1 




'S'31 


^32 


1 




- R33 


R11R33 


^33 - 




L '-'33 


-^33'S'33 


'-'33 J 



(3.3) 



We have that the Ihs of (|3.2p are meromorphic functions. It then follows that the matrix product 
on the rhs of these equations must also be meromorphic and must match any poles and residues 
found on the Ihs. Due to the poles in we can have double poles as well as single poles in ()3.2p . 
Prom the first column in the second equation in (|3.2p . at C = C2k^ obtain 

[R2ixt]{Ctk)=^ and 

Prom the second column and p.4p . assuming that S^^ does not vanish at (.^k^ obtain the single 
equation 

[xt + S32 4](Ctk) = 0^ (3-5) 
and from the third column, we obtain 

[W](C2+fe) = and 

Prom the above conditions, there are two different options in which these conditions may be satisfied. 
Thus we have a bifurcation where two different solutions are possible at each such common zero, 
which we now describe. 

3.2 UHP - Option UA 

Por the first solution, let us take the option to have 523 = = R21 at ^2^^, which we will call Option 
UA. Then the first equations in (|3.4p and (|3.6p are trivially satisfied while the second ones, along 
with (|3.5p . reduce to 

[^12 l^t + Xt] iCtk) = ' [^13/"^ - l^t] iCtk) = ' (3-7) 

where /i^ has a simple pole at (" = ^2^^, whose residue is X2'(C2^fc) / ^lif -^^'^ ^^^^ option, there is 
only one independent Jost solution at ^2^^- (For the adjoint Jost functions, there are two.) 



^+ - t;+7 {R'2ixt + R2ixt') + 



as 



-^11 '5'31 



f^3 



(3.4) 



-^13'S'33 ,,+ 



R^ 



+/ 



('S'23X^ + 'S'23X^') + l^t 



(C 



2±J 



(3.6) 
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3.3 UHP - Option UB 

For the other solution, the remaining option is to take 5*32 = = R12 at ^2^^- Then p.Sp and the 
first equations in (|3.4p and (|3.6p give that xt ™ust vanish at C^^- Thus 

X^(C2tfc) = 0, (3.8) 
while the second equations in (I3.4p and (13. 6p reduce to only one condition 



l^f^tiCtk) - -^f^tiCt,) + -^4iCt,k) = , (3.9) 

^ll,k '-'33, fc '-'33,A: 

where by I'Hopital's rule we have made use of iJi^iCtk) ~ xfiCtk) I ^iik- ^^^^ option, we have 
two linearly independent Jost solutions at ^2^^ (while the adjoint only has one). 

3.4 Bound state structure - LHP 

Now let us turn to the LHP and the structure of the Jost solutions when and R^^ have a 



common zero. Continuing as before but with the first equation in p.2p instead, from the first 
column we obtain 



[^21X21 (C2~fc)=0. 



Ml - ^ (•521X2 + ^21X2 ) + p-/ /^3 
-"33 -"-33 



(C2.) = 0, (3.10) 



from the second column and (|3.10p , assuming that Rfi does not vanish at a; ' obtain 

[X2 +^32A^3](C2:fc)=0, (3.11) 

and from the third column, we obtain 



^^/^r - ^ (^23x2- + ^23X2-') + M3- 

'-'11 -^11 



(C2;,) = 0. (3.12) 



[«23X21 iC2,k)=0, 

3.5 LHP - Option LA 

Continuing as before, under the option 5i2,fc = = i?32,fc at (2k^ have that X2 (C2~fc) ™ust vanish. 

From I'Hopital's rule it follows that A*2 (^2"^;) ~ ^2'(C2k) I ^33k S^'^^S /^2 ^° analytic at ^2"^- 
The other equations give only the one condition 

VMr(C2"fc) - (C2~fe) + ^3 (C2"fc) = . (3.13) 



'-'11, fc '-'11, fc -"'33,fc 

and there are two linearly independent Jost solutions for this option in the LHP (with the adjoint 
having only one such solution). 

3.6 LHP - Option LB 

Under the option 5*21, = = -R23,fc at C,^^.-, we find that there is only one independent solution for 
this option. The other two Jost solutions are given by 

(C2~fc) + -R32,fc H{^2,k) = ' Mr(C2"fc) - 'S'si.fc (C2"fc) = , (3.14) 

Here we have again that has a simple pole at C2k^ whose residue is X2^(C2^fc) / ^33' fc- 
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3.7 LDRs with Resonant Solitons 



Given the above, one may extend the LDRs given in Ref. [30] to include resonant sohtons. To do 
this, one would start with Eqs. (3.08)-(3.11) from that reference, which are still correct, even in 
the case of common zeros between the transmission coefficients. Starting from there, in order to 
reduce those Eqs. (3.08)-(3.11) to a minimal set of linear independent relations, it is necessary to 
remove all bound state <I>-type Jost solutions in favor of the ^-type Jost solutions. As was done 
for non-common zeros in Ref. [30], one can extend those equations by using the above options 
to eliminate the bound state $-type Jost solutions. The only change to the previous form of the 
LDRs, will be additional sums over any common zeros. (As we shall point out later, by treating the 
common zero case as a limit of different zeros approaching each other in Eqs. (3.15)-(3.19) of Ref. 
|30j . and by including the four possible options, one can obtain the same results as given below.) 

A priori, in the common zero case and without any symmetry on Q, one must allow all four 
options to coexist. To express this, we need to extend the notation in Ref. [30]. We shall use 
^2a k common zeros for option UA and LA with k = 1,2,..., A'^^ where is the number 

of these common zeros. Similarly ^ will designate the common zeros of options UB and LB 
where k = 1,2, . . . , with being the number of these common zeros. Meanwhile the current 
notation used in ()2.14p - ()2.19p will continue to be used for those zeros which are not common with 
any in the other transmission coefficient. 

Using the above options, it is straight forward to give these additional sums. One obtains for 

9(C) > 0, 



1 

2TTi 



nC'-C 



s 



11 



R 



33 



s 



11 



k=i (Cii,fc ~ C) 5*11 



^ ^''j^e2iCii,k)E2i{Cn,k) 



J2 ~\ 



S2I k R3I k 

-l-®2iC2a,k)E2liC2a,k) " ®3 (C2a,fc)^31 (C2a,fe) 

it. 



S- 



ll.k 



33,fc 



^2 

E 



^^©3 (C26,fc)-^3l(C2i,,fc) 



1 (^26,^: 'S'li^fc 



dC 



2vri JtzC'-C 



etEi2{C')--^Q3Es2{C' 



R 



33 



fc=l ''ll,A; -"-11, fc fc=i y^Z'i.k ^) -"-33, fc 



(3.15) 



^7^63 (C33,fc)^32(C33,fc) 



N. 



N- 



+ E 



T <^ti^tk)El2{Cta,k) + E 77-- 7-^ ^ 03-(C3-3,.)i?32(C2-,,.) , (3.16) 



1 C2a,k -( ^ 



11, fc 



'33,fc 
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and for $j(C) < 0, 



©2^(0 



1 



dC 



2vri JnC'-C 



it 



R 



■33 



E 



iCii,fc)^12(Cii,fcj +2^77^ 



1 Cii,A,. C Riij 



7l (^33,k -^33,fc 



7^ ©3 (C33,fc)^32(C33,fc) 



1 R 



fc=l ''2a,fc ^11, 



^2. 



1 i? 



N+ 

"33 



+ 



2vri Jn C'-C 



Sl3. 
•^33 



5'23 



33 



^7^03(C3-3,fc)i?32(C2-,,,), (3.17) 
'33,A; 



Rl2, 



^11 



1 S 1 5' 

XI 77+ TTT 02'(C33,fc)^23(C^3,fc) " X 77+ TTT 0+^ '^l(C2a,k)El3{Q 

j^^l \^33,k ^> '-'33, fc k=l ^^'2a,k ^) '^33,k 



2a,k) 



AO 

k=l ^'^2b,k ~ 



^^'^ Qt iC2b,k)El3iC2hk) - -^P^®2iC2b,k)E23{C2b,k) 



„+/ \'^2b,k'^^^\'^2b,k' q+f 

-"■11, fc '-'33, fc 



(3.18) 



These equations will be referred to as the "LDRs with resonances" or "extended LDRs" while those 
of [30] will be referred to as the "regular LDRs" or just plain "LDRs". 



3.8 Case of Non-compact Support 

The same relations as above could also be obtained in the case of non-compact support. We shall 
briefly outline here the approach that one would use. Starting from the relation (2.31) in Ref. [30j 
and the definition of the x states (the Fundamental Analytical Solutions, or FAS) in (2.15) - (2.19) 
of the same reference, one can obtain x"*" ■ X^"*" = Diagonal [Rfi, Rii S^^, 5-^] where is the 
adjoint FAS. (A similar expression also exists in the LHP.) Now at a non-common zero of either 
R^^ or 5^3 in the UHP, one has that there will be two linearly independent components of x , 
as described in [30]. However, if we have a common zero of these two transmission coefficients, 
then we have that x"*" ' X^~^ becomes a zero matrix, as well as the middle diagonal component 
becoming a double zero. From these conditions, one can deduce that at a common zero, must 
have either only one or only two linearly independent columns. (At the same time, the adjoint 
solutions, x^^7 correspondingly must have only two or only one linearly independent solutions.) 
These two conditions can then be seen to give rise to the above two options. The only difference 
would be that the coefficients in ()3.7p - (13. 9p will become arbitrary coefficients, generally unrelated 
to the components of R and S (since the off-diagonal components of these matrices do not generally 
exist off the real axis, in the non-compact support case). Once these options are established based 
on the required linear dependences and order of the zeros, then the same LDRs would result. 
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4 The Resonant (1, l)-Soliton Solution from the Extended LDRs 



Under the anti-hermitian symmetry, one can obtain a (1, l)-soliton solution when there is only one 
common zero between Ru and 533. First, take the scattering data to be in the form of Option A. 
Here we have = 1 = iV^ and whence k = 1 only. Thus the sum can be omitted. Furthermore 
the k subscript can also be deleted since it is no longer needed. Then Eqs. (j3.15p - (|3.18p become, 
for '=s{C) > 0, 



0^(0 



1 



and for 9(C) < 0, 



{C2a-0 

etiC) = 

02(0 = 



[C2102 (C2a)i^2l(C2;) " ^31 0^ (C2"a)^31 (C2"J] 



+ 



c 



12 



(2a c 



+ 



c 



12 



0^ 



C2a C 
Ci3 



0^(C2+)^12(C2+), 



C2a C 



Qti(2a)En{C 



2a J 



Let us define 



^12= ^^12(C2+a) 



521 

7^3l 



C21 
C31 



^2l(C2"J 
^3l(C2"a) 



gni2 gjC2a('^i-'^2)(a;-a;i2) 



V* 

g-nsi g-jC^„(Ji-J3)(a;-x3i) 



513= ^^13(C2+) 



7^31, 



From the above, we find the solution of (|4.ip - (|4.4p to be 

2?? 



0^(0 
0^(0 



+ 



03" (C) = 

where 77 is the imaginary part of (2a ^^'^ 



+ 



+ 



iC2a - C)Da 

2??7^l2 



■i{D - 1) 
—521 
^31 
1 

-iS2i 
^7^3l _ 

-1 
iS2i 



I)a = l + |7^12|' + |7^3l|' 



(4.1) 
(4.2) 

(4.3) 
(4.4) 

(4.5) 
(4.6) 
(4.7) 
(4.8) 



(4.9) 
(4.10) 
(4.11) 



(4.12) 
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From (j2.12p we then obtain the components of the potential matrix 



Qi2 = i{Ji - J2)2ri^ , (4.13) 
Qi3 = -i{Ji - h)^!!^ , (4.14) 
Q23= -ii J2 - J3)2v^^ , (4.15) 



with the other three components following from the symmetry relation on the potential matrix. 
Again, we have only two positions and phases to freely specify. The position and phase of Q23 can 
be considered to be driven by the other two. 

In the case of Option B we get essentially the same structure but with some indices interchanged. 
Taking 

7^l3= ^Ei3{Ctb)= e^T« e<^.(^i-^3)(x-xi3) ^ (4I6) 
531= ^^3i(C26) = n,, (4.17) 
7^32= ^Es2{C2b)= g-<,-,(J2-J3)(--32) ^ (4.18) 



523= ^E23{C+) = ^32, (4.19) 



one obtains for the potential components 



Qi2= i( Ji - J2)2j7^i|^ , (4.20) 



Qi3 = i{Ji - J3)2r/^ , (4.21) 
Q23 = -i{J2 - J3)2v^ , (4.22) 



where 

Db = 1 + lUrsf + . (4.23) 

Again we have one component being driven by the other two. In this case, it is Q12 which can be 
considered to be driven. 

The asymptotics of these equations are much easier to handle than those of (|2.22|) - (|2.29|) . 
In Option A, for X31 << X12, we have two NLS solitons present, with the one in Q23 far to the 
left and the one in Q12 far to the right, with Q13 vanishingly small. As we reverse this and take 
X12 « X31, then we find that the original two NLS solitons have vanished and a single NLS soliton 
has appeared in Q13. 

For Option B, if we take X13 << X32, then we have two NLS solitons present, with the one in 
Q12 far to the left and the one in Q23 far to the right. As we reverse this and take X32 << 2:13, then 
we find that the original two NLS solitons have vanished and a single NLS soliton has appeared 
in Qi3- From these asymptotics, one sees clearly that these two options are distinctly different 
solutions since the spatial orientation of Q12 and Q23 become reversed. 

There is one case where the solutions of these two options do overlap. This is when TZi2 = = 
^23, giving that Q12 = = Q23 with only Q13 remaining nonzero, for each option. This case is also 
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the third case wherein a single sl{2) sohton can occur in this system. As was mentioned earher, 
the (0, l)-soliton and the (1, 0)-sohton are NLS sohtons, with one in Q12 and the other with one in 
Q23- Each one of these NLS sohtons only require one zero in one transmission coefficient. Here we 
see that Q13 can also contain a single NLS soliton, but to obtain this sl{2) soliton, one has to insist 
that there be a common zero in two different transmission coefficients, a more difficult objective to 
achieve practically. 

Prom the above, one can also argue the instability of the resonant soliton. The NLS solitons 
in the (0, l)-soliton and the (1, 0)-soliton solutions are known to be stable. Small variations in any 
part of their scattering data gives small variations to the NLS solitons. On the other hand, the 
existence of the Qis NLS soliton is predicated on more than one coefficient being exactly zero. 
Let either one of these coefficients become shifted, then one has the case where the two common 
eigenvalues are no longer exactly equal. As discussed in Section 2.4 and as seen from (j2.22p - (j2.29p . 
the asymptotics will be distinctly different from those found above for the resonant soliton. 

5 The Resonate Soliton as a Limit of One Eigenvalue Approaching 
the Other 

Here we will briefly discuss how the resonate soliton can be obtained as a limit of a zero of R^^ (C) 
and a zero of 5';^(C) approaching each another. As we have already seen, one has two choices as 
to how this can be done. Once that is realized, then one must choose one of two options and 
it becomes straight forward to proceed in that manner and obtain the extended LDRs from the 
regular LDRs. We will leave that as an exercise to the reader. 

The same thing can be done with the general (l,l)-soliton solution, ()2.26p - (]2.29p . Looking at its 
structure, one sees if one chooses C12 to be proportional to {Ct ~Ct) (|2.22p and also renormalizes 
X12 to take this into account, then as the eigenvalues approach each other, the solution degenerates 
into the Option A case. Similarly by taking C23 to be proportional to {C,f — C,'^) in (I2.25|) . and 
renormalizing X32, then one obtains the solution for Option B. 

6 Conclusion 

There were two reasons for this study of the soliton solutions of the 3x3 operator. The first was to 
obtain a better understanding of the resonant soliton and the scattering data which gives rise to it. 
The second was to gain a better understanding of the possible soliton solutions in the general n x n 
eigenvalue problem. As has been shown here, the resonant soliton can be obtained from either the 
LDRs or the (M, iV)-soliton solution, by taking the limit of two eigenvalues in different transmission 
coefficients approaching each other. The general (M, A^)-soliton solution for the 3x3 operator has 
been constructed in Appendix A from the LRDs. The approach taken there was one that could be 
extended to the n x n case, upon being given the LRDs for the n x n case. (For the case of the 
nxn operator, one would include more algebraic equations of the form given in ()A.3P -( El4l) . which 
could then be solved in the same manner outlined in Appendix A.) Of course, as n increases |49j . 
the number of transmission coefficients will also increase, resulting in even more complex forms of 
resonant soliton solutions, as well as increased options for the number of independent Jost functions 
at any given resonant eigenvalue. 

The simpler soliton solutions of the 3x3 eigenvalue problem have been presented, including 
the resonant soliton. In the general case, there are four different transmission coefficients, two in 
the UHP and two in the LHP. The zeros of these give rise to general sl{3) solitons and generally 
will express themselves as multiple sl{2) solitons in the appropriate components of the potential 
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matrix Q. However whenever any of these sohtons are near to one in the other wave, the corner 
components of Q could also be excited. 

There are several symmetries which could be applied to the 3x3 potential matrix and each 
one would generally have different soliton structures |24[ |47] . For the anti-hermitian symmetry, the 
eigenvalues and normalization coefficients in the LHP become equal to the complex conjugates of 
their counterparts in the UHP. In this case, one can designate the ,3/(3) soliton solutions by (M, A^) 
where M and N are the number of zeros in the two transmission coefficients in the UHP. The 
simple case of a (1,0)- or (0, l)-soliton solution is identical to an NLS soliton. These solitons will 
only be found in one of the two off-diagonal components which are directly next to the diagonal 
of Q. Moving up to the next case of a (1, l)-soliton solution, as long as the two transmission 
coefficients do not have eigenvalues close to one another and as long as the NLS solitons do not 
have their positions close together, the (1, l)-soliton solution is essentially a sum of two (1,0)- and 
(0, l)-solitons. 

When the magnitude of the difference between the two eigenvalues become significantly smaller 
than their imaginary parts, the s/(3) (1, l)-soliton approaches a (nonlinear) resonance. In this 
realm, the s^(3) soliton is found to have a presence in all off-diagonal components of Q. As the 
positions of the two NLS solitons come closer together, the amplitude of the s/(3) soliton in the 
(5i3 wave increases while the components in the other waves Q12 and Q23 correspondingly decrease. 
However the component of the soliton in the Q13 wave never fully forms into a NLS soliton. As the 
difference between the eigenvalues becomes smaller and smaller, the Qis component comes closer 
and closer to forming a complete NLS soliton. This is the regime of a near nonlinear resonance 
wherein the soliton in Qis never does fully form. 

When the eigenvalues become exactly equal, one has the fully nonlinear resonant case. Here 
the (1, l)-soliton solutions are best given by (I4.12p - (l4.15p and (I4.20p - ()4.23p . However as pointed 
out earlier, they can also be obtained as a limit process from (j2.22p - (j2.29p . The range of solutions 
here vary from widely separated NLS solitons in Q12 and (523 with essentially no component of the 
s/(3) soliton in Qiz, to the NLS solitons in Q12 and (523 asymptotically being completely absorbed 
into a fully formed NLS soliton in (5i3- The asymptotic limit of this latter solution provides us 
with the third sl{2) soliton. This asymptotic limit for ()4.12p - (j4.15p follows upon taking C12 = 
and the same can be obtained from (j4.20p - (|4.23p upon taking C23 = 0. 

As shown earlier, there are two different structures for resonant solitons, depending on whether 
there is one (Option A) or two (Option B) independent Jost solutions at the eigenvalue. The 
general solution for Option A of the resonant (1, l)-soliton has only one independent Jost solution 
at the eigenvalue and is given by ()4.12p - ()4.15p . Equations ()4.20p - ()4.23p give the solution for 
Option B which has two independent Jost solutions at the eigenvalue. As pointed out earlier, the 
solutions for these two options are distinctly different except for the one common solution where the 
normalization coefficients C12 and C32 both vanish. When that happens, the solution degenerates 
into the form of a single NLS soliton in Q13 (and Q^i) with all other components of Q zero. Here 
we have a single NLS 1-soliton solution resulting from two zeros, each in a different transmission 
coefficient. As mentioned earlier, in this resonant case, viewed from the evolution of the 3WRI, the 
solution for one option is just the time-reversed solution of the other. 

The reason this resonance is a "nonlinear" resonance is that it requires an exact equality between 
two eigenvalues, each of which is a zero of a different transmission coefficient, and each of which 
is independent of the other. From the direct scattering problem, the positions of these zeros 
are nonlinear functions of the potential matrix Q. From the inverse scattering problem, it is the 
eigenvalue which determines the spatial shape of the amplitude and phase of a soliton solution. (The 
normalization coefficient determines only the position and overall phase of the soliton solution.) 
When these two eigenvalues are equal, then the shape of the solitons in Q12 and (^23 are such that, 
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in the language of the 3WRI, as these two waves coUide, they will coherently and constructively 
interfere with each other, transferring all their energy into the third wave, Q13. Recall that for 
plane waves in the 3WRI, the resonance conditions are that the sums of two of the wave-vectors 
and the sums of two of the frequencies must equal those of the third. However for plane waves, one 
never achieves a complete conversion since the inverse processes will limit and prevent complete 
conversion. Thus the resonant soliton is a solution to the 3WRI problem of how to shape and phase 
the envelopes of two colliding beams, Q12 and Q23, such that one could achieve a total conversion 
into the third wave, Qis- The condition for this nonlinear resonance and complete conversion is 
the equality of the eigenvalues. Similarly, any fully resonant (every eigenvalue for one wave being 
paired with another in the other wave) {N, A^)-soliton solution would also do the same. 

A Appendix A: Derivation of the general soliton solution; non- 
common zeros 

We obtain the (M, A^)-soliton solution in the absence of any symmetry. We assume no common 
zeros between the M zeros and the zeros and that each M zero and each N zero is simple. One 
can substitute the expressions for the 02 's from (I2.15P and (I2.17P into the remaining expressions 
([2T3]1 and ([2^8]) to get 
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where the J-''s and ^'s are defined in the obvious manner. Now, using these expressions and 
evaluating Qf{C) at C = Cn_fc, 
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for A; = 1, 2, ... , A''^^, and evaluating 63 {() at C = ^33 ^, 
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A; = 1, 2, ... , A^33, which together are a system of N^^ + ^'33 vector equations for the 0]'"(C^ ^)'s 
and ©3 (C33fc)'s- Solving this system, and plugging them back into the original expressions (|A.3P 
and ()A.4p . we recover the closed form solutions for 0^(C) and 63 (C), as well as Qf{C,). Explicitly, 
we may write the system in matrix form, to wit: 
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where the +l's occur with the diagonal terms only. Then, (assuming Ai is invertible) the general 
solutions for 0^((^^ ^) and 03'(C33A:) O^y Cramer's rule) 
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where A4k is the matrix formed by replacing the A;th column of M with a and Mj^^+^f^ is the 

matrix formed by replacing the (A'^xi + k)tli column of A4 with a. Arranging all terms properly, we 
obtain 
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Under this framework, these results scale up for the soliton solutions of the general n x n operator 

EH. 
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